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Abstract 

The behavior of scalar perturbations on superhorizon scales during the reheating stage is 
investigated by replacing the rapidly oscillating inflaton field by a perfect fluid obtained 
by spacetime averaging and the WKB approximation. The influence of the energy transfer 
from the inflaton to radiation on the evolution of the Bardeen parameter is examined for 
realistic reheating processes. It is shown that the entropy perturbation generated by the 
energy transfer is negligibly small, and therefore the Bardeen parameter is conserved in a 
good accuracy during reheating. This justifies the conventional prescription relating the 
amplitudes of quantum fluctuations during inflation and those of adiabatic perturbations 
at horizon crossing in the post-Friedmann stage. 
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§1 Introduction 

In the inflationary paradigm cosmological large scale structures such as galaxies and their 
distribution are thought to be formed from seed density perturbations produced by quan- 
tum fluctuations of an inflaton field]!], H ]. In this scenario one can in principle determine 
the statistical properties of the present large scale structure of the universe by calculating 
the amplitude and the spectrum of the seed perturbations and tracing their evolution, if 
the fundamental laws of nature are specified. The former task is relatively easy and the 
result can be put into a simple formula which is valid for a wide variety of the inflaton po- 
tential. For simplified evolutionary universe models, it is also the case for the second task 
as far as the linear evolutionary stage of perturbations is concerned. This simplification is 
brought about by the fact that a gauge invariant variable called the Bardeen parameter is 
conserved with a good accuracy for growing modes in Friedmann stages and the inflation- 
ary stage 0, H |, H ]. For example, if the reheating at the end of the inflationary stage is 
instantaneous, one can determine the amplitudes of perturbations when they reenter the 
horizon in the post-Friedmann stage from those of the quantum fluctuations on the Hubble 
horizon scales during inflation simply by matching the values of the Bardeen parameter. 

This powerful conservation law of the Bardeen parameter still holds for smooth transi- 
tion from the inflationary to the Friedmann stage, provided that the equation of state of 
the cosmic matter changes slowly with the cosmic expansion and that entropy perturba- 
tions can be neglected@ ]. In realistic models of reheating, however, it is not clear whether 
these conditions are satisfied or not because matters with different dynamical properties 
such as the inflaton field and radiation coexist possibly for a long while. Further, according 
to the perturbation theory of multi-component systems, the energy transfer processes, for 
example, from the inflaton to radiation, themselves may produce additional entropy pertur- 
bations if the energy transfer rate depends on the energy densities or the cosmic expansion 
rate. For example, it is proposed recently that a parametric resonance of the coherent in- 
flaton field and other massless scalar fields may contribute as the dominant energy transfer 
process in the early phase of reheating]^, [7], ||, [5]]. In this case large entropy perturbations 
can be produced by the perturbations of the energy transfer rate since this rate for the 
parametric resonance is very sensitive to the amplitude of the inflaton oscillation and the 
cosmic expansion rate. 

In the present paper we investigate this problem in detail and examine whether the 
Bardeen parameter is well conserved or not during the reheating phase by evaluating the 
entropy perturbations produced by realistic reheating processes with the help of the gauge 
invariant formalism for mult i- component systems^ ]. On the basis of the result in our 



previous workJT0| ] that the behavior of perturbations in the stage dominated by an oscil- 
latory inflaton <fi coincides with that of a perfect fluid obtained by a spacetime averaging 
and the WKB approximation except for a sequence of negligibly short intervals around the 
zero points of 0, we consider a system consisting of gravity, radiation and a perfect fluid 
corresponding to the inflaton. 

The paper is organized as follows. First in the next section we give the basic assump- 
tions and explain their motivations. In particular the behavior of a scalar field in the flat 
spacetime which decays to a massless scalar particles via the standard one particle decay 
process(Born decay) is briefly analyzed to prove the preservation of coherence of the decay- 
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ing inflaton field. This analysis is also used to determine the structure of the energy transfer 
term. In §3 the evolution equations for perturbations during reheating are given by apply- 
ing the gauge-invariant formalism for the perturbations of a multi-component system to 
our system. Then in §4 the amplitudes of entropy perturbations produced during reheating 
are estimated with the helps of these evolution equations to show that the Bardeen pa- 
rameter is conserved with a good accuracy for perturbations corresponding to the present 
large scale structures. §5 is devoted to conclusion and discussion. The gauge-invariant 
perturbation theory of a multi-component system used in §3 is reviewed in Appendix A. 
The definitions and formulae are given in a generic form in order to correct errors in the 
corresponding equations in the original article ]. In Appendix B a general method to get 
upper bounds on solutions to a first-order differential equation system is explained. It is 
used to find upper bounds on the entropy perturbations in §4. 

Throughout this paper, the natural units c = h = 1 are adopted and 8nG is denoted as 
k 2 . Further the notations for the perturbation variables adopted in the article [|| ] are used 
and their definitions are sometimes omitted except for those newly defined in this paper. 



§2 Fundamental assumptions and preliminary consid- 
erations 

In order to investigate the evolution of perturbations during the reheating phase, we must 
treat a coupled system of the inflaton, matter and gravitational field. In this paper, as 
the cosmic matter, we only consider radiation whose dynamics is described by the energy- 
momentum tensor of a perfect fluid 

T(r)nv = (Pr + Pr)U(r),Ji(r) v + g^P r , (2.1) 
Pr = \pr, (2.2) 

where and in the followings the tilde denotes perturbed quantities. Further we assume 
that the inflaton is well described by a classical and coherent real scalar field minimally 
coupled with gravity, and its energy-momentum tensor is given by 

f Wflu = d^d u 4> - \g„» [(V0) 2 + 17(0)] . (2.3) 

Due to the decay of the inflaton these energy-momentum tensors are not conserved sepa- 
rately: 

= Own = -Qv> ( 2 - 4 ) 

^y T (r)fi = Q{r)n = Qfi- (2-5) 



2.1 Coherence of the inflaton 

One subtle point in the assumptions above is the requirement of coherence on the inflaton 
field. If the interaction of the inflaton field and radiation is neglected, no problem arises 
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when we assume that it is described by a coherent classical field. On the other hand, 
when the interaction is taken into account, the consistency of this assumption is no longer 
obvious because the inflaton field by itself is a quantum field and its decay into radiation 
is a quantum process in general. Further such quantum nature should be properly taken 
into account when one determines the structure of the energy transfer term Q^. 

Though we cannot justify this assumption for a generic case, it seems reasonable at least 
in the case in which the potential of the scalar field is quadratic in from the following 
observation. 

Let us consider a massive scalar field on a flat background interacting with a massless 
scalar field x which plays the role of radiation in the realistic situations. For simplicity we 
assume that the interaction of these fields is cubic and their Lagrangian density is given 
by 



£ = — 
2 



(dcf)) 2 + m 2 (j) 2 + (dx) 2 + 2/x0x 2 • (2-6) 



Let us quantize these scalar fields by the standard canonical quantization and introduce the 
creation and annihilation operators, (Ap,Ap) and (B^B^), in the Schrodinger picture 
by 

4>{x) = E -7=2= (A p + Alp) e*P x , (2.7) 
p ^J2Vuj p 

= -i E Op - a -p) ( 2 - 8 ) 

X (^'(A k + A^ t (2.9) 



= -*EV^( A fe + A U) e<fc "* ( 2 - 10 ) 

where and n x are the conjugate momentums for and Xi respectively, and 



u p := ^\p\ 2 + m 2 . (2.11) 

Then these canonical variables are represented on the Fock space Ti = TL<f, <8> H x spanned 
by the basis vectors 

|n, m >= -^(AT) % ® f E ^=r(^) m ) (2-12) 

where Q = ® Q x is the Fock vacuum defined by 

ApVL^ = 0, B k tt x = 0. (2.13) 

In this representation the free part of the Hamiltonian for each component is written 
in the standard diagonal form: 

:H^.= J2^ApA p , :H x :=J2\k\BlB k . (2.14) 
P k 
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Hence Q is the ground state of the free part of the Hamiltonian, but not of the total 
Hamiltonian. Though this implies that Q is unstable against the time evolution, it is a 
convenient cyclic vector for constructing coherent states. 

Let us consider a state $(£) which is represented at the initial time t = as 

$(O)=$ ®ft x . (2.15) 

Then by the perturbation theory, taking account of the mass renormalization, the expec- 
tation value of the annihilation operator for is calculated as 

< <S>(t)\A p \<S>(t) >= e"i < A p (t) > +0(/i 3 ) , (2.16) 

where Ap(t) is the operator defined by 

A p (t) := Ape^ 1 —Alp cos u p t, (2.17) 

UlpTp 

< X > denotes the expectation value of X in the state $0, and r p is the life of 0-particles 
with a momentum p given by 

T * = lifc < 2 - 18 » 

The expectation value of the creation operator is just given by the complex conjugate of 
the above expression. 

Next we consider the expectation values of products of creation and annihilation oper- 
ators. It is convenient to define << X » t by 

« X » t :=< $(t)\X\$(t) > - < n(t)\X\Q(t) >, (2.19) 

where Q(t) is the solution to the Schrodinger equation with the initial condition Q(0) = Q. 
This subtraction of the expectation with respect to Q(t) is to eliminate the effect of the 
instability of the perturbative vacuum. In this notation, the expectation values of products 
of creation and annihilation operators are given by 

« ApA q » t = e -* ( i + i } < : A p {t)A q {t) :> +0(V) , (2.20) 
« ApAq » t = e -* ( i + i } < : Ap(t)A q (t) :> +0(^ 3 ) . (2.21) 
In particular if we take the initial state $^ as a coherent state given by 

= exp Yl l c p| 2 j ex P (l2 c v A P^j n <f» ( 2 - 22 ) 

the expectation values of Ap and Ap, and their products with respect to this initial state 
are simply given by 

< A p > = c p , < Ap > = c p , (2.23) 

< ApAq > = c p c q , < ApAq > = c p c q . (2.24) 



Hence from Eq. (pl6Q , Eq.(^20|) and Eg .(g2|) we obtain 



« X(zOny) »*-< >< >, (2.25) 

up to 0(/x 2 ) where X and K are any of <j) and 7T^. 

This result shows that the coherence of the scalar field is preserved by its Born decay 
at least within a time shorter than the decay life r. Though the above argument does not 
apply to the energy transfer by the parametric resonance, it is reasonable to assume the 
coherence of the inflaton in that case as well because the parametric resonance occurs only 
when the inflaton field has a good coherence. 



2.2 Energy- momentum transfer term 

The simplified model analysis in the previous subsection can be used to determine the 
structure of the energy-momentum transfer term as well. 

In that model the divergence of the energy-momentum tensor of the quantum 0-field is 
given by 

where the second term on the right-hand side is the mass counter term. Calculation of the 
expectation value with respect to $(t) in the order /z 2 yields 

/r <: d^(t, x) ]T — L= K p {tyP- x :> , (2.27) 



duTfa. = fi X 2 d^ - Sm 2 ^, (2.26) 



16vr p UpW 



0(t, x) := £ 1 (A p e~ luJ ^ + Alpe^) e l P x , (2.28) 



where 



it p (t) := -%^ v (Ape-^ 1 - A_pe^ 1 ) . (2.29) 

In the case in which the initial coherent state $^ contains only particles with small mo- 
mentums, as in the case of inflaton, this expression is approximately written as 

Q(c/>)n = -7 <: d^(t, x)n(t, x) :> 

~ -7^ < 0(t, x) > < it(t, x) >o= -7^0vr, (2.30) 

where 

7 :=-^-. (2.31) 
167rm 

This result suggests that in curved spacetimes in general the energy transfer term for 
the Born decay has the form 

= -Q Wli = d^ffa n v d„4>), (2.32) 
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where n M is some timelike unit vector which coincides with the unit normal to the constant 
time hypersurfaces in the spatially homogeneous case. Though we cannot determine this 
unit vector, its choice has no effect in the framework of the linear perturbation theory 
for the following reason. In the linear perturbation, since n J is a first-order quantity, the 
perturbation of n^d^cf) depends only on Sn° as 

Sipfdjt) = d 5(f) + 5n°d <j). (2.33) 

On the other hand, for the same reason, 5n° is determined only by the perturbation of goo 

as 

Sn o = _1 ^po_ (2 34) 

2 #00 

Hence the freedom of n M has no effect. 

The above argument on the structure of Q M may not apply to the energy transfer term 
for other processes such as the parametric resonance. However, since <9 M is the unique 
vector field constructed from <p in the inflaton dominated stage, it is reasonable to assume 
that has the same structure as that given in Eq.( |2.32 ) for such cases as well, though V 
may depend on higher derivatives of (f) 



2.3 WKB approximation and replacement of the scalar field by 
a perfect fluid 

From the discussion so far we can formulate the evolution equation for perturbations dur- 
ing the reheating stage by applying the gauge-invariant formalism for a general multi- 
component system to the system consisting of the classical scalar field, radiation and the 
gravitational field. However, the equations obtained by this procedure is rather difficult to 
analyze because some of the terms in the equations becomes very large periodically when 
vanishes. 

In order to avoid this difficulty and make the problem tractable, we utilize the result 
of our previous paper [|T0| ]. It is shown there that during the stage in which the rapidly 
oscillating classical scalar field dominates the energy density of the universe and its behavior 
is well described by the WKB form 

4> = F(S,x), (2.35) 

where S is a rapidly oscillating phase function, the behavior of superhorizon scale pertur- 
bations coincides with that for a perfect fluid system obtained by a spacetime averaging 
of the energy- momentum tensor of Tjt^ v over Hubble horizon scales except for a sequence 

of negligibly short intervals around the zero points of 0. This implies that we can replace 
the rapidly oscillating scalar field by a perfect fluid in the investigation of the behavior of 
superhorizon perturbations as far as the perturbation variables averaged over the oscilla- 
tion period are concerned. On the basis of this result we consider a perfect fluid instead of 
treating the classical scalar field directly. 
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The energy-momentum tensor of the perfect fluid corresponding to the classical scalar 
field is given byfllOl ] 

ffa v = (p f + Pf)u1 f) u if> + PSS, (2.36) 
P f = Wfpf, w f = (2.37) 

where we have assumed that the potential of the scalar field is given by a simple power-law 
function 

U = -\<j)\ n . (2.38) 
n 

We assume this power law form throughout this paper, pf and uf^ is represented in terms 
of the original scalar field as 

Pf = ^ < (V0) 2 >- ^ < (dsF) 2 > (V < S >) 2 , (2.39) 

u {) = 5 < ^ > 1 ; (2.40) 
im [(V<5>) 2 ]^' V ; 

where < X > represents a spacetime average of X over the Hubble horizon scales. This 
approximation is good if the parameter defined by 

e 2 ■= < ( VF ) 2 > (2 41) 

is much smaller than unity. Since e represents the ratio of the cosmic expansion rate H 
to the oscillation frequency of the scalar field, this condition is satisfied in the rapidly 
oscillating phase. 

In order to formulate the perturbation equations for this perfect fluid and radiation, we 
must rewrite the spacetime average of the energy-momentum transfer term Q2.32j ) in terms 
of the fluid variables. If we write < > as for simplicity, it must have the form 

Qfi = < S > R(< (f> > 2 , n u d u <S>). (2.42) 

From the argument in the previous subsection we can take < S > as h^. Hence, from 
the relativistic virial theorem [|i~0| ] 

< u$) >= _I < (V0) 2 > [1 + 0(e)] , (2.43) 
n 

it is simply written as 

Q, = u (f) ,Q, Q = G(p f ). (2.44) 

For example, for the Born decay in the model considered in the previous subsection, G is 
simply given by 

2n 

G(Pf) = ^tfrpf ( 2 - 45 ) 
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On the other hand, for the energy transfer by the parametric resonance, Q should be 
modified as 

Q = G{p f ,±V lt u 1 }), (2.46) 

where V^Wj/3 represents the expansion rate of the <fi = const hypersurface, and coincides 
with the Hubble parameter H in the unperturbed background. This dependence arises 
because the duration of the parametric resonance for each mode of massless fields coupled 
with (j) depends on the cosmic expansion rate. Though V^Wy is of order e in the WKB 
approximation scheme, it may not be neglected because Q has a strong dependence on it 
for the parametric resonance decay. 



|3 Evolution equations for perturbations 



Under the assumptions given in the previous section, we can easily write down the evolution 
equation of perturbations during reheating by applying the gauge invariant perturbation 
theory of a multicomponent system to the current system. Basically a perturbation of this 
system is described by the gauge invariant density contrasts and the shears of the 4- velocity 
for the inflaton fluid and radiation. But in order to investigate the dynamical behavior 
of the Bardeen parameter, it is more convenient to choose variables which respect the 
decomposition of perturbations into the adiabatic modes and the entropy modes. Hence, 
as the basic variables, we adopt the curvature perturbation the total shear velocity 
V, and the quantities representing the difference of the density contrasts and the shear 
velocities between the inflaton and radiation, Y and Z, defined by 



Y 



PfPr 
P 2 



s 



fri 



PfPr dH 



P A 



k 



Vi 



fr- 



(3.1) 



Note that Y and Z become zero at the beginning and at the end of the reheating stage 
when the energy density of radiation or the inflaton is negligible. 

With the help of the general formulae given in Appendix A, we can easily write down 
the evolution equations for the basic variables by specializing the space dimension to 3. 
First, from Eq.( |A.16]) , ( |A.17|) , (|A.31|) , and (|A.33|) , the evolution equations for $ and T are 
written in terms of these variables as 



-I K (l + w)T, 




T 



1 1 ( , l i 



»f- 3 Y, 



(3.2) 



(3.3) 



where T> = a(d/da) = d/d(loga). On the other hand, by taking account of Eqs. ( |A.21| ), 
(|A.22p , and (|A.28[) , Eq.( |A.37|) and (|A.38|) reduce to the following evolution equations for 
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Y, Z: 



VY + 3[ 1 - + w f -2w)Y = -(±-) Z- 3 -®-±±^E cf 
\3 ' / \aH) AHpl + w s c} 



VZ + 



29. _L_ h + i Wf ) 9jl (A.) », 

3Hpl + w f V 4 } ) I K \aH J 

3 /i 

' a + w)I K - + 3(w f - 2w) + (1 - $w f )-£ 



a 2 H 2 





+ Hh 


1 




~ h 








+§( 





2 4 



1 



Pr 



3 3 1 + Wf pf J 
Y 



Z 



r + -(l+W f 

3) I K \aH . 



In particular, subtracting Eq. 



from Eq. 



p 2 1 + w 
we obtain 



(3.4) 



(3.5) 



{ ' a 2 H 2 3\aH 



4 1 + w f 
"3(1 + w) 2 



i i 



i 



I K l + w h{ W ' h ' + 3 h 'l t 



(3.6) 



where $ — T denoted as ( in the article [ITT] ] is referred to as the Bardeen parameter 
from now on[fJ ]. Taking into account that $ — T and $ are of the same order, and 
that the spatial curvature K is practically zero because of the inflationary expansion, we 
immediately confirm from this equation that the Bardeen parameter $ — T is conserved on 
superhorizon scales, if the entropy perturbation Y is negligibly small. In the next section, 
we will investigate in how good accuracy the Bardeen parameter is conserved by evaluating 
the amplitude of the entropy perturbation generated in realistic reheating processes. 

In the evolution equations of the entropy perturbation, the perturbation of the energy 
transfer rate works as the source term. When the energy transfer rate Q depends only on 
Pf, the perturbation of the energy transfer rate is expressed in terms of the basic variables 
as 



E, 



cf 



A 



G p f (pf)pf 

G(Pf) 



G Pf(pf)pf^ + w f 
G{p s ) 1 + w 




3 P/ 



(3.7) 



where G n , denotes the partial derivative of G with respect to pf. If the energy transfer 



'pf 



rate depends also on the cosmic expansion rate V nU f as 



Q 



t 1 f 

G{~ Pf ^V^f) 



(31 
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the following term should be added to the right hand side of Eq. i 



G H (p f ,H)H fl k 
G(p f ,H) \3aH f 



K 



2 R 2 



G H ( Pf ,H) H 
G(p f ,H) 



k 



3 \aH 



T 



T 



4 1 



k 



P 



K 



9 1 + w \ aH J pj 



a 2 H 2 



T 



(3.9) 



On the other hand, if V^w^, or V M w^ are adopted as the local Hubble constant, the terms 
to be added are given by 



G H {p f ,H) H (I k y K T N 



G(p f ,H) \3aH 

G H ( Pf ,H) H 
G{p f ,H) 



a: 



> H 2 



i fjr 

3 [aH 



T 



K 



2 R 2 



T 



(3.10) 



and 



Vr 



K 



G H (p f ,H) H fl k 

G(p f ,H) \3aH" a 2 H 2 



■T 



G H (p f ,H) H 
G(p f ,H) 



1/feV l l + w f 
3\aH 3 1 + w 



aH 



Pr 



K 



a 2 H 2 



T 



(3.11) 



respectively. 

In all of these expressions for E c f, all the terms in proportion to $ or T are multiplied 
by coefficients of order (k/aH) 2 . In the next section, we will show that this suppression 
factor makes the contribution of the entropy perturbations negligibly small, even if we take 
into account a possible dynamical growth of Y and Z. 



»4 Conservation of the Bardeen parameter during the 
reheating stage 



From Eq. ( |3.6| ), it follows that the entropy perturbation does not affect the conservation of 
the Bardeen parameter if the condition 



a <= da, 



\y\ < m 



(4.1) 



is satisfied, where a s and a e are the values of the scale factor at the start and at the end of 
the reheating stage. In this section we show that this condition are satisfied in the realistic 
chaotic inflation scenario whose dominant reheating processes are the parametric resonance 
and/or the Born decay. 

First we define the index g p and gu by 



9 Pf 



G p f Pf 
G ' 



9h 



GrH 
G ' 



(4.2) 



10 



Though g p . is bounded by unit from below, 



g P} > i, (4.3) 

for realistic reheating processes, it may become much larger than unity in the stage in 
which the parametric resonance is effective. It is difficult to evaluate the upper bound 
on g p . in this stage because we have poor knowledge on the strong parametric resonance. 
However, since the analysis in the weak parametric resonance || ] shows that G behaves as 
G ~ exp(0(<7 P/ )) when g p ^ 1, it is expected that g p does not exceed unity by many 
orders of magnitude. For example, according to the recent numerical investigation taking 
account of rescattering of produced particles|T2| ], the effective value of g Pf does not exceed 
100. In contrast to g p , gu does not have a definite sign and may becomes negative. 
However, from the analysis of weak parametric resonance, it is expected that its absolute 
value \qh\ is & t most of the same order as g Pf . Hence, when the distinction of these is not 
important, we use 

g := max{g Pf ,\g H \}. (4.4) 

Note that as the Born decay dominates in the energy transfer, gu vanishes and g P} = g 
approaches unity. 

In order to evaluate the amplitude of entropy perturbations, it is convenient to decom- 
pose the reheating stage into the following four substages: 

1) a s < a < a\ 

2) cii < a < a 2 



G 

H Pf 




G 

-n—g > i 

Hpf 


> G 


Hp,' 


1 > G 

1 ~ rr 9, 

Hpf 


1 ~ J. 


G 

IT- >>L 
Hpf 





G 



HPf 

4) a 3 < a < a e > 1. (4.5) 

Hpf 

The first substage corresponds to the period during which an explosive energy transfer 
occurs by the parametric resonance. As the amplitude of the inflaton oscillation decreases 
and the parametric resonance gets less effective, G/Hpj becomes less than unity. If g is 
much greater than unity in this phase, the second substage appears. On the other hand, 
if the Born decay already dominates at that time and g — 1, this stage does not appear 
and the system goes directly to the third substage, during which the Born decay is the 
main process of reheating but it is still slower than the cosmic expansion. Finally as the 
cosmic expansion rate decreases with time, G/Hpf becomes greater than unity again and 
the reheating completes. This is the last substage. 

We evaluate the upper bound of the amplitude of the entropy perturbation in each 
stage by utilizing the technique explained in Appendix B. We assume K — henceforth. 
First we put Eqs. ( |3.4[ ) and ( ^.5[ ) into the matrix form 

VX = CtX + 5, (4.6) 
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where 2-column vectors X and S, and 2x2 matrix O are defined by 
( Y \ 



X 



ft 



^ — 1 — 3wj + 6w 

~ H Pf 9pf 



) 2 ("1- 



:9h- 



i 



3Hp f ^ n l+wf' 
3 i 9„„ g„„ . ( i 0„,,.\h r 



ff P/ ll" 7 / ^ 3> h ^3 ft J / 

/ -lG_ n 4- 1 G 4+3w f \ 
2 _ffpyp/ T 6 Hp 1+Wf \ 



S :-- 



2 l — 3Wf pfpr 

' 9 1+w p 2 
1 G „ 1+w 



+ 



4 Hp9H 1+Wf 




(4.7) 

In order to estimate the upper bound of := fl + f2', we decompose it into a sum of 
three matrices as 



ft 



HI 



-2 - Qw f + 12w 








-3 + 9w - Qw f - (2 - Qw f )^ J ' 



f2 



_H"2 



1 + 



G 



3# P/ 9H\ +Wf 

U 
h 



i 



) \ 



if3 



(off) + 3Hp f 9H 1+Wf ) 







-^{(2 W/ + + 
Then the maximum eigenvalue of these matrices are given by 

A m i = max | — 2 — 6u>/ + 12w, —3 + 9w — 6wf — (2 — Qwj)-^- 



(4.8) 



A 



m2 



G 1 \ 4 p r 1. .py 



\ m3 = - min <i 2 



'if,/'" 

respectively. Hence the maximum eigenvalue of f2# is bounded by the sum of them, 
A m : = A m i + A m2 + A m3 



G 



2w f + t)^ + ^ 



1 p r 



3/h 31 + Wf h j ' 



(4.9) 
(4.10) 
(4.11) 
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< N- 



A 

aH 



1 + 



G 



-9h- 



f G G 
mm < 2— — g p 



SHpf 1 + Wf J 

4\ p f 8 1 



w/ h 



(4.12) 



Here A" is the maximum value of the sum of A m i and the two last terms on the right hand 
side of (|4 . 1 Q|) . Since this sum becomes maximum at p r = or pf = for fixed Wf, N is 
given by 

5 

N := max{2 - 5w f ,&w f - -}. (4.13) 



The source term is evaluated as 




l$l 



(4.14) 



taking account of the fact that $ is of the same order as T. 

First in the first stage, since A m is negative and its absolute value~ G/Hpf is much 
larger than unity, Eqs.( p.l"0|) and ( p.!4j) in the Appendix B yields 



\X(a)\\<0 {9 (± 



l$l 



(4.15) 



Next in the second stage, A m is bounded as A m < since (k/aH) 2 g is practically much 
smaller than unity. Hence, by using Eq.( B.l(j|) in the Appendix B and by taking account 
of the order of magnitude of \\X\\ at the end of the previous stage, ||X(ai)||, we obtain 



i*(«)ii < 



.1 



N/2 



O gi 



k 



l$l 



,ai / \ \ aH 

where gi is the value of g at a — Oj. In the third stage, in the same way we obtain 



(4.16) 



\X(a)\\ < 



a \ 
aiJ 



N/2 



Oft 



A 

aH 



l$l 



(4.17) 



Finally in the fourth stage, since the eigenvalue of O is negative, and its absolute value 
increases exponentially in time, the influence of the previous stage on X(a), corresponding 
to the first-term in Eq.( [B.10|) , is rapidly erased, and ||X(a)|| settles down to a value 
determined locally by ||5i|| as 



\X(a)\\ < 




(4.18) 



as seen from Eqs.( ETTU ) and ( B.14 ). This estimate is actually a quite weak one. The actual 
value of ||X|| decreases exponentially in time in this stage due to the suppression factor 
Pf/p in the definition ( |3.1|) . Hence Y and Z take nonnegligible values only in stages when 
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the inflaton and radiation coexists, and vanish as soon as the energy transfer completes 
effectively. 

These estimates on the upper bound of ||X|| can be used to obtain stronger bounds on 
Y and Z. To see this, let us write the second row of Eq.( f4.6|) as 

VZ = Q Z Z + S z , (4.19) 

Q z := -| + |«- 3^ - (1 - 3w f )^ (4.20) 

S z :=0(Y ) + 0^)\y (4.21) 

and regard this as the evolution equation for Z. Then since fl z is bounded as 

n z <N z -o(J^J, (4.22) 

N z := max j-3 W/ , ^ - 1)} , (4.23) 

using the values of the upper bound obtained in the previous analysis on \Y\(< \\X\\) 
and applying Eqs.( |B.l"0|) and (p,14| ) to the above equation on Z, we obtain the following 



stronger upper bound on \Y\ in the first and the fourth stage: 



Hp f k 



2 



We can apply the same method to the evolution equation for Y obtained from the first 
row of Eq. ( |4.6j ) , 

VY = VL Y Y + S Y , (4.25) 

n Y :=-l-tou f + Gw- -=—g Pr (4.26) 

Now Qy is bounded as 

n Y <N Y --£-g Pf , (4.28) 

N Y :=\l-3w } \. (4.29) 

Hence we obtain in the first stage 

m«)l<0| [^] 9], (4.30) 
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assuming that 



in the second stage 



\Y(a)\ <0[max<l,S. [£) m (AV J ] I __ I 4 , (4.32) 




and in the third stage 



2 



Putting these estimates together, we finally obtain 

— \Y\ = / — |Y| 

a A NY n( Ji ^ 
O max < 1, oi — 




< [f J O^l^fJ ( ^| Hl^l 1*1- («4) 

Taking account of N/2 > Ny when n > 2, we can conclude that the entropy perturbation 
Y does not affect the conservation of the Bardeen parameter if 

is satisfied. Since g\ is a monotonic function of pf and it is expected that it does not exceed 
unity much, we drop it henceforth. 

The condition ([4.35 ) gives a lower bound on the energy density p(a e ) at the end of 



reheating. Let us evaluate its order of magnitude for realistic values of the physical param- 
eters. First we introduce the parameter y defined by p(a e ) 1//4 = KFGeV. This parameter 
is related to the value of k/aH at the end of reheating for perturbations corresponding to 
the present cosmological structures of 10 Mpc scales by 

(A) A 10— . (4.36, 

Here we have assumed a(t eq ) / a(t ) ~ 10~ 4 where t eq is the equality time. In order to gener- 
ate density perturbations consistent with the observed anisotropy of the cosmic microwave 
background ST/T ~ 10~ 5 , the energy density at the time when the relevant perturbations 
cross the Hubble horizon in the inflationary stage should be given by 



,V4 _ io 16 GeV, (4.37) 
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which is approximately equal to p(a s ). Hence from the inequality 




1/[3(1- 



mm < Wf 



(4.38) 



we can see that the condition ( f4.35|) holds, if 



y > < 



576wf — 356 

36wf 
221w f 



13 
19 



5w f - 10 



for 



for 



Wf 



w f 



> 



< 



3' 



(4.39) 



is satisfied. 

This condition on y is a strict one from a mathematical point of view. However, it 
seems to be too strong from a physical point of view because Y takes nonnegligible value 
only around x := p r /pf ~ 1, while the above estimate is based on the upper bounds on A m 
and fly that correspond to the values at x = or x = oo. Hence a practical condition is 
obtained by replacing N and Ny in Eq. (|4.35|) by the values of A m and fly at x — 1. It is 
expressed as 



y > 



lOOwf + 404 

29w/ + 57 
135u$ + 370w f + 299 

18w^ + 65w f + 43 



for 



for 



Wf 



Wf 



> 



< 



3 ■ 



3 ' 



(4.40) 



These conditions are depicted in Figure 1. In this figure the solid line and the dotted 
line correspond to Q4.39 ) and ( 4.40 ), respectively, and the Bardeen parameter is well con- 
served in the region above these lines. This figure shows that reheating does not affect the 
conservation of the Bardeen parameter if it terminates before the primordial nucleosyn- 
thesis. Hence we can conclude that in all the realistic models based on chaotic inflation, 
the Bardeen parameter stays constant in a good accuracy during reheating. This jus- 
tifies the conventional prescription relating the amplitude of adiabatic perturbations at 
horizon crossing in the post-Friedman stage to the value of the Bardeen parameter in the 
inflationary stage. 



§5 Discussion 

In this paper we have investigated the evolution of perturbations during reheating taking 
account of the effect of the energy transfer from the inflaton to radiation, by replacing 
the inflaton field by a perfect fluid obtained by the spacetime averaging and the WKB 
approximation. By evaluating the amplitudes of entropy perturbations generated during 
reheating and their influence on the adiabatic component of perturbations, we have shown 
that the Bardeen parameter is well conserved during the reheating stage as well as in the 
inflationary stage and the post-Friedmann stage for realistic models. Though we have 
considered only the parametric resonance and the Born decay as the dominant reheating 
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Figure 1: The energy density at the end of reheating (p 1 ^ 4 
inflaton potential n 



KFGeV) vs power of the 



processes for dehniteness, the conclusion holds rather generally since the arguments are 
insensitive to the details of the models. 

Of course our analysis does not exhaust all the possible models of inflation. In particular 
we have only considered the case in which the inflaton is described by a single component 
field. Though a simple multi-component extension does not seem to change the conclusion 
as far as the fluid replacement of the inflaton fields during the reheating stage gives a good 
approximation, some subtlety may occur if scalar fields with very tiny masses coexist with 
the inflaton fields. For example, if there exists a scalar field which affects physical param- 
eters controlling the reheating processes such as particle masses or coupling constants of 
physical particles but has no dynamical effect by itself during reheating, it may produce 
large entropy perturbations and affect the Bardeen parameter. In particular, as our anal- 
ysis suggests, this possibility may become important if such a field affects the parameters 
controlling the parametric resonance processes. In such multi-component systems, entropy 
perturbations produced during reheating may survive after reheating and have important 
effects on the present universe, as discussed by Yokoyama et al. as a mechanism to gen- 
erate the baryon isocurvature perturbation 1 13 ]. These problems in the multi-component 
extension are under investigation. 
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Appendices 

A Gauge invariant Perturbation Theory of Multi Com- 
ponent Systems 

In this appendix, we recapitulate the basic definitions and equations of the gauge invariant 
perturbation theory for multicomponent systems developed in the article|5| ]. The purpose 
is twofold. The one is to explain the notations and to provide the basic equations used in 
the text. The other is to correct errors of the corresponding equations in |5] ]. 

The Einstein equations and the equations of motion of a perturbed multi- component 
system are given by 

& v = k 2 ^^ (A.l) 

V ^ T (a) M = Q(a)ft, (A. 2) 

where Q( a )^ represents the energy-momentum transfer term for a component (a), which 
satisfies 

E <9(a)„ = 0- ( A -3) 
(a) 

In a spatially homogeneous background with metric 

ds 2 = -dt 2 + a 2 dal, (A.4) 

where da 2 denotes a constant curvature space of dimension n with a sectional curvature 
K, these equations reduce to 



rji {A 2 2k 2 K 

H =(-) = ~( V\P 2' 


(A.5) 


T>p = —nh, 


(A.6) 


T>p a = -n(l - q a )h a , 


(A.7) 


d 

V:=a—, 
da 


(A.8) 



Here 



h = p+P and h a = p a +P a denote the entalpies of the total system and the (a)-component, 
respectively, and q a is defined by 

q a = -^-- = (Q a ,0,...,0). (A.9) 

nh n ri 



For scalar perturbations, the linear perturbation equations for Eq. ( A. 1|) are written in 



terms of gauge-invariant quantities A, V, $, T, and II representing the perturbation 
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amplitudes of density, velocity, curvature, gravitational potential, entropy and anisotropic 
stress, respectively, as 



k 

VA - nwA = -C K (1 + w) — V -in- l)C K wU, 

aH 

VV + V = 

k T k 

* H 

aH aH 



-A + r 



l + w 1 + w 



n ~ l k w 
n aH l + w 



nix a?H 2 

a 2 

(n - 2)$ + * = -k 2 — pll 

A; 2 



(A.10) 

(A.ll) 
(A.12) 
(A.13) 



Here w = Pj p, c 2 = Pj p and 



Ik:=1 + 



K 



a 2 H 2 



C K := 1 



nK 



(A. 14) 



In applications it is often more convenient to rewrite these evolution equations in terms of 
$ and T defined by 

T := ?fv, (A.15) 



as 



X?$ + (n - 2)$ = -^I K (1 + w)T - ^-^(^) 2 /^n, 
n K 

vr + -i K (i + w )r-—r 



+< 



_( n _ 2 ) + ^^-( — 

1 + w n aH Ik 

n(n — 1) aH 

2 /, 



n-l„ 
a <^a 



l + w 
l + w 



n. 



(A.16) 



(A.17) 



On the other hand, from Eq.( [A.2| ), the perturbation equations for each component are 
given by 



V[p a A a ] + np a A a 



-C 



A 



k 



K — V a + (n - 1)(1 - q a )PJl 
aH 



aH. 



-nq a [P a T a + c a p a A a ] + nh a q a E a + n{\ - q a )—h a F a 



VV a + V a = ^^> + ' 



c 



aff aff \ 1 + w r 



-A r 



1 + W a 



n — 1 k „ w a 



A' " 



n aif 1 + w a 



-n Q + f q . 



(A.18) 



(A.19) 
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Here E a and F a are the gauge-invariant quantities representing the energy and the mo- 
mentum transfer rate for the component (a) at its own rest frame. These quantities are 
related to the corresponding quantities at the rest frame of the total system, E ca and F ca , 
by 



F n 



E a + ^(V a -V), 



F a + nq a (V a -V). 

These quantities must satisfy the following constraint obtained from Eq. (|A.3| ): 

^ t QaE C a 0, ^ h a F ca 0. 



(A.20) 
(A.21) 

(A.22) 



Eqs.( [A~T0D , flATH ), flOSp , (|A~T9|) and determine the dynamical behavior of 

the multicomponent system completely. However it is often more convenient to use the 
variables representing the relative perturbations among components defined by 



S a /3 : 

v aP -■ 
n a/ 3 



■P 



1 + W a 1+Wp' 

V a - V p , 



l + w a 



Wp 



1 + Wp 

wp 

l + Wp 



n/j, 



(A.23) 

(A.24) 
(A.25) 

(A.26) 



l + w a 

q a E ca — qpE c p, (A. 27) 

F ca - F c p. (A.28) 

The dynamical variables for each component are written in terms of these relative variables 
and those for the total system as 



E a p 

Fa/3 



1 + w a 
V n = 



l + w + 2 ^ h^' 



w. 



1 + W a 
W n 



r a 
n, 



w 



1+w 
w 



T - T re i) + "7^ r "7> 



(A.29) 
(A.30) 
(A.31) 
(A.32) 



1 + w a 1 + w 

where r re i is the contribution to T from the relative entropy perturbations and given by 



PT 



rel 



EE 



h a hp 
h 



(l-q a )(l-qp)(c 2 a -c 2 p) x 
A c/3 



(l + w a )(l-q a ) (l + wp)(l-q p )_ 

2 



^EE 

Z a p 



(A.33) 
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The evolution equations for these relative variables are derived as follows. First, by 
subtracting Eq. flA.lll ) from Eq. (|A.19|) , we obtain 



A 



k 



v(Va - v) + (K - v) = _(4 - 4)— + ^fcis, 



cry 



k W 



aH 1 + w ' c ' ' aif 4^ h L 
+F a + n4(l-g a )(F a -F). 

Next from the relation 



n — I 



n 



A, 



l + w a 1 + w a ' l + w a p a k 
Eq- dOD , and Eg flOl , it follows that 



(A.34) 



(A.35) 



1 + w 



A 



+ to 



+n(l - g Q ) - T - nq a (l + c 2 a ) ^ "r 1 ^ 



-n- 



1 + w 

— r Q - (n - l)C K (l - q a 



h 

w 



-II + nq a E c 



1 + w« 1 + w 

From these equations we obtain the following evolution equations for S a p and V a p\ 



(A.36) 



n , 



VS af3 + -[q a (l + cQ + qp{l + 4)}S a p 

fi h 

+f + 4) - + 4)] E f + sw 



it 



a/7 



- nT aB + nE, 



af3 



-n[(q a - qp)c s + q a {l + cj - ^(1 + c p )]- 
-w(?a - 9/9)7- — r + (n - l)Cjf(g a - g^) — — n 

1 + W 1 + U7 



+ U> 
If 



(A.37) 



n , 



VV aP + V a/3 --[c 2 a + cj- q a (l + 4) - g^l + 4)]^ 



77 



r 2 2 
2 i C a _ C f3 



aH [ ° a Cf3) l + w 1 2aH 



q a (l + c 2 a ) + q p (l + 4)] \ (V<n + V Pj ) 

A + o^f(4 + 4)s a/3 + 2^(4 - 4) E + ^7) 



afl" 



a/3 



n — 1 



(A.38) 
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B Upper bound on the growth rate of solutions to a 
first-order differential equation system 

In this appendix we explain a general technique to evaluate an upper bound on the norm 
of the solution X to the first differential equation system 

jX = n{t)x + s{t). (b.i) 

Here X and S are n column vectors, and ft is an n x n matrix. If we define the norm of 
the solution X by 

\\X\\ 2 :=X j X, (B.2) 



it obeys the equation 



^(\\X\\ 2 ) = X^n H X + S^X + X^S, (B.3) 
at 



where fin is an hermitian matrix defined by 

:= fi f + Ct. (B.4) 

The least upper bound of the right-hand side of this equation is determined by the 
maximum eigenvalue of fin- However, in many cases this simple method is not practical 
because it gives a complicated expression. In such cases it is useful to decompose into 
a sum of hermitian matrices ft Hi, (1 < i < k) as 

k 

Q H = Y,^m, (B.5) 

i=i 

so that the maximum eigenvalue of each Qhi has a simple expression. Since the hermitian 
matrix can be diagonalized by some unitary matrix and its eigenvalues are real, we 
obtain 

X^Q m X < X mi X^X, (B.6) 

where X m i is the largest eigenvalue of the hermitian matrix ft Hi- Hence by applying the 
Cauchy-Schwartz inequality 

TiA ] B < (TrA t A) 1 / 2 (TrB t S) 1 / 2 (B.7) 

to Eq.( [B.3| ), we obtain 

|(||X||)-^||X||-||5||<0, (B.8) 

where 

k 



A m := ^2 ^rni- (B-9) 
i=l 

By integrating the inequality ( B.8|) , we obtain 



X 



< eKpUZfldt?) \\X(0)\\ + j[Vexp [jl^Y^dA \\S{t')\\. (B.10) 
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Notice that the right-hand side of the inequality (p.lOQ is a monotonically increasing func- 
tional of A m (s), (0 < s < t). Therefore, even if we replace A m by a larger value, the 
inequality ( |B.10|) still holds. 

In the case X m (t) < and |A m (t)| ^> 1, we can derive a simpler estimate on ||-X"||. To 
see this, let us introduce x{t') defined by 

x(t') := - f \ m (t")dt" . (B.ll) 
Jt> 

Then the second term in the right-hand side of Eq. (|B.10| ) is written as 

jjVexp (U^m(t")dt") \\S(t')\\ = £dxexp (-^ (B.12) 

where 

A= f\\ m (t')\dt'. (B.13) 
Jo 

Hence, if ||S||/|A m | is a slowly varying function of t and A > 1, we obtain 
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